Observation of Periodic Orbits on Curved Two— dimensional Geometries 
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We measure elastomechanical spectra for a family of thin shells. We show that these spectra can 
be described by a "semiclassical" trace formula comprising periodic orbits on geodesies, with the 
periods of these orbits consistent with those extracted from experiment. The influence of periodic 
orbits on spectra in the case of two-dimensional curved geometries is thereby demonstrated, where 
the parameter corresponding to Planck's constant in quantum systems involves the wave number 
and the curvature radius. We use these findings to explain the marked clustering of levels when the 
shell is hemispherical. 
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When studying spectra of quantum systems, semiclas- 
sics, in particular periodic orbit theory, provides a pow- 
erful connection to the dynamics of the analogous clas- 
sical system [IHl]. Considerable progress was recently 
made in understanding generic statistical properties of 
quantum chaotic systems 5J. Furthermore, spectra of 
specific systems can be constructed from periodic orbits, 
prominent examples are the Hydrogen atom in a strong 
magnetic field 6], the Helium atom [7] as well as regular 
and chaotic two-dimensional billiards [3] . Microwave 
experiments in flat cavities are most useful because, in 
two dimensions, the Helmholtz equation for the electrical 
field formally coincides with the Schrodinger equation. 

Do the concepts of quantum chaos carry over to other 
wave phenomena? — Elasticity, i.e. mechanical vibra- 
tions, is a particularly interesting testing ground. The 
governing equations and the boundary conditions are 
different from the ones in quantum mechanics. More- 
over, different modes (pressure and shear) are present and 
propagate with different velocities, in case of anisotropy 
even depending on the directions. Hence, a transfer 
of quantum chaos ideas to elastic systems is a worth- 
while endeavor in its own right. The bulk waves in 
three-dimensional aluminum blocks have the same sta- 
tistical features as known for quantum chaotic sys- 
tems Remarkably, even much more subtle fea- 
tures have been measured and understood in a frame- 
work transferred from quantum chaos: parametric statis- 
tics |10j , transport and localization properties in three- 
dimensional blocks [TT], the statistics of elastic dis- 
placements, i.e. "wave functions" , in two-dimensional 
plates [12] and Wannier-Stark ladders in quasi-one- 
dimensional elastic systems [13]. 

In this contribution, we measure spectra of thin elastic 
shells and explain certain characteristics by periodic or- 
bits. This is, to the best of our knowledge, the first exper- 
imental identification of periodic orbit features on curved 
shells. We have two goals. First, we show that periodic 



orbit theory can be applied to the shells. This is non- 
trivial because of the curvature and because the modes 
are a combination of flexural and in-plane fields. They 
are in general described by a system of partial differen- 
tial equations of high order jl4j , very different from the 
Schrodinger equation which is the starting point in quan- 
tum chaos. Second, we will use this insight to explain a 
striking clustering effect in the spectra. These findings 
are not only of conceptual but also of practical impor- 
tance. Metallic shells are ubiquitous in technology and 
every day life, ranging from auto bodies to micro-electro 
mechanical systems. Since the calculation of spectral fea- 
tures from the wave equation is often tremendously com- 
plicated, an understanding in terms of simpler geometric 
quantities, the periodic orbits, might be of considerable 
interest. 

The shells we employ are objects of revolution as shown 
in Fig. [l] They form a family of constant mid-surface 
area parameterized by the opening angle Oq measured 
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FIG. 1: (Color online) Shells of revolution with opening angle 
6o and curvature radius R. 

from the axis of revolution, whereas R is the 9q depen- 
dent curvature. For 9^ — 7r/2, we have a hemisphere, 
the angle 6*0 = formally corresponds to the planar disk. 
Figure [2] displays the aluminum shells used in the exper- 
iment, the opening angles are 6*0 — 0°,25°,52° and 90°. 
Each shell is 2 mm thick although for the hemisphere 
three thicknesses were studied. The disk has a diameter 
of 80 mm, the other shells were made in such a way that 
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FIG. 2; (Color online) Aluminum shells used in the experi- 
ment. From left to right, the opening angle is 60 — 0° , 25°, 52° 
and 90°. 



all mid-surface areas are equal. The shells were carved 
from a solid aluminum block to avoid any internal strain 
which inevitably occurs when producing curved objects 
by bending a plate. It was carefully made sure that the 
cut defining the boundary was along the radius, as seen 
in Fig. [1] This is important to prevent coupling between 
flexural and membrane-like modes. 

We used the experimental setup developed in our pre- 
vious studies [TOl [12] . It has a very high resolution, the 
quality factor Q = //A/ is typically about 10^ where / 
and A/ are position and width in frequency of a given 
resonance, respectively. We accumulated data in the fre- 
quency range kHz < / < 800 kHz. In the measurement, 
the shells are only supported by three 1 mm ruby spheres 
which minimize the elastic coupling to the rest of the 
setup. The acoustic coupling to the air is significantly re- 
duced by putting the setup into a vacuum chamber whose 
inside is held at 10"'^ Torr. Flexural (bending) modes 
and membrane-like (stretching and shearing) modes are 
excited. The resulting spectra are displayed in Fig. |3] 
A visual inspection immediately reveals that the spectra 
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FIG. 3: (Color online) Measured spectra for the shells with 
opening angles 9o = 0° , 25°, 52° and 90° from top to bottom. 

for the disk and the shells with opening angle Oq < tt/2 
look rather similar. For the hemisphere with = ^/2, 



however, the levels are structured in clusters which are 
almost equally spaced. The clustering was already found 
in Ref. |l5j by numerically solving the equations for thin 
shells [13]. In the following, we will use periodic orbits 
to give a clear and intuitive explanation of this striking 
effect. 

Our reasoning will be in the spirit of semiclassical anal- 
ysis for quantum systems. We demonstrate that the 
in-plane excitations can be described by a sum over pe- 
riodic orbits, i.e. by a trace formula akin to the ones in 
quantum chaos [lH5- This goes well beyond previous 
work on the celebrated whispering gallery modes [2], on 
the ray-description of seismic waves |16j as well as on 
the identification of periodic orbits in elastic spectra of 
three dimensional systems [101 ET] . Trace formulae are 





FIG. 4: (Color online) Geodesic pentagram and diameters on 
the spherical cap. 

different for regular and chaotic systems. If the wave 
equation describing shell vibrations were scalar, all our 
shells would be integrable. The complexity of the wave 
equation for thin shells [H] modifies that picture. The 
wave equation can be separated into an angular part de- 
pending on the angle of revolution and a "radial" part 
depending on the angle 9 measured from the axis of rev- 
olution with < < Oq. In the "semiclassical" approxi- 
mation [18| the waves are approximated by the motion of 
a fictitious particle. The main insight is that this motion 
takes place on the geodesies of the shells as in Fig. |4] The 
role of Planck's constant compared to a typical action is 
here played by 1/kR, where k is the wave number. To 
leading order, the fiexural and the membrane-like mo- 
tion decouple, i.e. yield two separate equations of mo- 
tion for the fictitious particle. Hence, the fiexural motion 
is integrable. There are two degrees of freedom on the 
shells and two constants of motion, the energy and the 
angular momentum with respect to the axis of revolution. 
The membrane-like motion, however, has two polariza- 
tions, one longitudinal (L, pressure) and one transverse 
(T, shear) mode. They are always coupled upon reflec- 
tions at the boundary. In this sense, the membrane-like 
motion is not integrable. 

The measured spectrum in Fig. [3] lowest plot shows 
evidence of eigenvalue clusters which are separated by 
a slowly increasing spacing in frequency. The first clus- 
ters appear only after a gap. At medium to high fre- 
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quencies the spacing appears constant and for thin hemi- 
spheres agrees weU with the spacing coming from purely 
T-polarized orbits. Furthermore, the extracted discrete 




FIG. 5: (Color online) Time spectrum of spectral fluctua- 
tions containing Longitudinal, Transverse (a) and Flexural 
segments of orbits (b). 

spectrum show fluctuations well described by certain 
time periods, see Fig. [Sj corresponding to fixed spacings 
in frequency which we attribute to periodic orbits in the 
following. From shell theory, we therefore focus on the 
in-plane deformation field for which dispersion is linear 
and results in orbit actions linear in frequency. We shall 
also touch briefly on the flexural mode as its dispersion 
approaches linear at the higher frequencies probed in our 
experiment. 

The in-plane theory results from ignoring the flexural 
motion in the Kirchoff-Love shell theory [HKTS] and leads 
to 
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when dropping curvature terms in the ray limit. The 
constants E and v are the elastic modulus of extension 
and the ratio of Poisson corresponding to the plane stress 
approximation used in plates |14j . Up to further curva- 
ture terms, this curved version of Navier-Cauchy's elastic 
equation is rewritten by re-expressing the covariant vec- 
tor laplacean VaV" in terms of two-dimensional curls 
and a gradient on a divergence |19] . Consequently, the 
elastic field can be decomposed in a gradient field and 
a curl field corresponding to longitudinal and transverse 
polarization each satisfying a curved scalar Hclmholtz 
equation: 

(A + z2)uW^0 (2) 



with i = L,T a polarization index in the following and 
Zj = kjR dimensionless wave numbers obeying zt = kz^ 
with K = cl/ct the ratio of propagation speeds in plates. 
On the sphere, deformations are therefore expanded 
using Legendre functions {P, Qjjf (cos 0) exp(im(/)) with 
n — Ij a polarization dependent angular momentum 
obeying z'j — + 1). We impose free boundary con- 
ditions: the integrated stress tensor across the thickness 
of the shell vanishes at a boundary. Thus, the normal 
components of the stress resultant vanishes: Ng^ = at 
9 = 71/2 with a = 9,4> [H]. The spectrum is then found 
by a two-by-two determinantal condition. 

As in the scalar problem [TB] we use the method of 
scattering quantization [20] and consider the dynamics of 
a ray of polarization j in the angular momentum variable 
Ij evolving under the condition that is constant. By 
elementary calculations we find this conservation law to 
be equivalent to the laws of reflection and refraction for 
the ray. For the derivation of a trace formula we therefore 
sum over m = = II sin 9 with 9 the incidence angle 
with respect to the normal. Except for the increase in 
complexity, we find a very similar result as in |18j : there 
is a scattering matrix composed of a propagation part 
and a reflection/refraction part. 

The propagation part over a single great arc of the 
hemisphere is found to be exp(— zttZ-,) with Ij « Zj + 1/2 
with the 1/2 related to the caustic phase shift. The free 
propagation over the sphere is then independent of m 
and no saddle point integration is needed as in the case 
of an opening angle different from 7r/2. 

The total m-dependence resides only in the reflection 
coefficients aq used for constructing an orbit q. For large 
I each reflection coefficient is asymptotic to the classical 
reflection coefficient and a slowly varying function of to, 
so the sum over to is well approximated by an integral 
and yields 

amaq — 



tr an 
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d9 cos9aq{9) = II Aq 



The final result of our calculations is that the oscillat- 
ing number of states dN in a frequency interval df coming 
from orbits is given to leading order as 

dN = pip{f)dj K.lL^Afp cos{-Krlp)dlp , (4) 



where is a closed orbit containing an r'th repeat of a 
prime ray sequence p and Ip = ulIl + utIt with nL^nr 
integers. We find it useful to group orbits with identical 
actions such as (LT)'^ and LLTT and denote the family 
L?T^ . In practice A for each family is calculated numeri- 
cally. Notably, fluctuations grow linearly with frequency 
to leading order. Rayleigh edge orbits are omitted as 
their contribution is of order /° only. 
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For completeness we include transverse rays with inci- 
dence angles beyond the critical angle for conversion for 
which the reflection coefficient is of unit modulus. Con- 
sequently, this part of phase space for T-polarized orbits 
leads to clusters of closely grouped states [21] ■ Ref . [H] 
states that Q should hold in general; here for the first 
time we present an explicit study in a ray splitting case. 
We give a parallel treatment to the flexural orbits using 
their corresponding reflection coefficients and disper- 
sion relation Here we settle for one extracted from 
the numerical solution of a whole spherical shell due to 
its larger spectral range. 

The simulated spectral fluctuations of the hemisphere 
used in our experiment have a period spectrum depicted 
in Fig. [5] the (a) spectrum shows clear evidence of in- 
plane orbits. The inclusion of the flexural modes in (b) 
gives trains of peaks with the main peaks close to those of 
the transverse. This is because at higher frequencies the 
flexural dispersion approaches the linear with the flexural 
speed approaching the Rayleigh speed which is just below 
the transverse speed. Furthermore, the clusters in the 
transducer signal we checked to agree well with the peaks 
in the flexural fluctuations. 

Why is the spectrum for the hemisphere so different? 
This is illustrated in Fig. |4j For all opening angles, all 
periodic orbits are geodesic polygons. Fig. |4] (left) shows 
the pentagram as an example. As the opening angle in- 
creases, the surface enclosed by the polygons becomes 
larger. At 9o = 90°, this surface is the entire hemisphere. 
Hence all polygons degenerate to the orbit on the base 
of the hemisphere. The diameter orbit shown in Fig. [4] 
(right) is the exception: there is only one such orbit for 
6*0 < 90°, because it must go through the north pole. At 
6*0 = 90°, however, all geodesic lines connecting opposite 
points on the base are diameter orbits. These drastic 
changes in the periodic orbit structure reflect the equally 
drastic changes in the spectrum when the opening angle 
reaches 90°. Formula Q gives the precise mathematical 
connection. 

In conclusion, we measured high-resolution spectra of 
a family of open spherical shells and saw clear evidence 
of in-plane behavior. For this we developed a trace for- 
mula for the density of states which agree well with ex- 
periment. The amplitudes of the fluctuations were not 
entirely as in experiment, still the location of the pe- 
riods of the experimental fluctuations agreed well with 
theory. Hence, we have presented a flrst identification of 
periodic orbit structures on curved geometries. From the 
simple form of the trace formula we expect it to be easy 
to generalize and apply to other wave equations on the 
hemisphere. Thus, when applied to the flexural mode, 
also the clusters of the transducer signal are described. 
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